Whole genome non-invasive prenatal testing in prenatal screening algorithm: clinical experience from 12,700 pregnancies.

BACKGROUND: A fast adoption of a non-invasive prenatal testing (NIPT) in clinical practice is a global tendency last years. Firstly, in Russia according a new regulation it was possible to perform a widescale testing of pregnant women in chromosomal abnormality risk. The aim of the study-to assess efficiency of using NIPT as a second-line first trimester screening test in Moscow. METHODS: Based on the first trimester combined prenatal screening results 12,700 pregnant women were classified as a high-risk (cut-off ≥ 1:100) and an intermediate-risk (cut-off 1:101 - 1:2500) groups followed by whole genome NIPT. Women from high-risk group and those who had positive NIPT results from intermediate-risk group were considered for invasive prenatal diagnostic. RESULTS: 258 (2.0%) samples with positive NIPT results were detected including 126 cases of trisomy 21 (T21), 40 cases of T18, 12 cases of T13, 41 cases of sex chromosome aneuploidies (SCAs) and 39 cases of rare autosomal aneuploidies (RAAs) and significant copy number variations (CNVs). Statistically significant associations (p < 0.05) were revealed for fetal fraction (FF) and both for some patient's (body mass index and weight) and fetus's (sex and high risk of aneuploidies) characteristics. NIPT showed as a high sensitivity as specificity for common trisomies and SCAs with an overall false positive rate 0.3%. CONCLUSIONS: NIPT demonstrated high sensitivity and specificity. As a second-line screening test it has shown a high efficiency in detecting fetus chromosomal anomalies as well as it could potentially lower the number of invasive procedures in pregnant women.

Gene Loss, Pseudogenization in Plastomes of Genus Allium (Amaryllidaceae), and Putative Selection for Adaptation to Environmental Conditions.

Amaryllidaceae is a large family with more than 1,600 species, belonging to 75 genera. The largest genus-Allium-is vast, comprising about a thousand species. Allium species (as well as other members of the Amaryllidaceae) are widespread and diversified, they are adapted to a wide range of habitats from shady forests to open habitats like meadows, steppes, and deserts. The genes present in chloroplast genomes (plastomes) play fundamental roles for the photosynthetic plants. Plastome traits could thus be associated with geophysical abiotic characteristics of habitats. Most chloroplast genes are highly conserved and are used as phylogenetic markers for many families of vascular plants. Nevertheless, some studies revealed signatures of positive selection in chloroplast genes of many plant families including Amaryllidaceae. We have sequenced plastomes of the following nine Allium (tribe Allieae of Allioideae) species: A. zebdanense, A. moly, A. victorialis, A. macleanii, A. nutans, A. obliquum, A. schoenoprasum, A. pskemense, A. platyspathum, A. fistulosum, A. semenovii, and Nothoscordum bivalve (tribe Leucocoryneae of Allioideae). We compared our data with previously published plastomes and provided our interpretation of Allium plastome genes' annotations because we found some noteworthy inconsistencies with annotations previously reported. For Allium species we estimated the integral evolutionary rate, counted SNPs and indels per nucleotide position as well as compared pseudogenization events in species of three main phylogenetic lines of genus Allium to estimate whether they are potentially important for plant physiology or just follow the phylogenetic pattern. During examination of the 38 species of Allium and the 11 of other Amaryllidaceae species we found that rps16, rps2, infA, ccsA genes have lost their functionality multiple times in different species (regularly evolutionary events), while the pseudogenization of other genes was stochastic events. We found that the "normal" or "pseudo" state of rps16, rps2, infA, ccsA genes correlates well with the evolutionary line of genus the species belongs to. The positive selection in various NADH dehydrogenase (ndh) genes as well as in matK, accD, and some others were found. Taking into account known mechanisms of coping with excessive light by cyclic electron transport, we can hypothesize that adaptive evolution in genes, coding subunits of NADH-plastoquinone oxidoreductase could be driven by abiotic factors of alpine habitats, especially by intensive light and UV radiation.

Detection of S-HBsAg Mutations in Patients with Hematologic Malignancies.

Multiple studies of hepatitis B virus (HBV) genetic variability and its relationship with the disease pathogenesis are currently ongoing, stemming from growing evidence of the clinical significance of HBV mutations. It is becoming increasingly evident that patients with hematologic malignancies may be particularly prone to a higher frequency of such mutations. The present report is the first extensive study of the prevalence of escape mutations in S-HBsAg, performed using isolates from 59 patients from hospital hematology departments with diagnoses of leukemia (n = 32), lymphoma (n = 20), multiple myeloma (n = 3), and non-tumor blood diseases (n = 4). The isolates were serologically examined for the presence of HBV markers and sequenced using either next-generation sequencing (NGS) or Sanger sequencing. Occult hepatitis B was found in 5.1% of cases. Genetic analysis of the region corresponding to S-HBsAg demonstrated an exceptionally high mutation frequency in patients with leukemias (93.4%) and lymphomas (85.0%), along with the prominent mutation heterogeneity. Additionally, more than 15 mutations in one sample were found in patients with leukemias (6.3% of cases) and lymphomas (5.0% of cases). Most of the mutations were clinically significant. The study analyzes the mutation profile of HBV in different oncohematological diseases and the frequency of individual mutations. The data strongly suggest that the NGS method, capable of detecting minor populations of HBV mutations, provides a diagnostic advantage, lays the foundation for the development of screening methods, and allows for the study of the virological and pathogenetic aspects of hepatitis B.

Complete plastome sequencing of Allium paradoxum reveals unusual rearrangements and the loss of the ndh genes as compared to Allium ursinum and other onions.

In this work the complete chloroplast DNAs of Allium paradoxum and Allium ursinum, two edible species of Allium subg. Amerallium (the first lineage), were sequenced, assembled, annotated, and compared with complete Allium plastomes of the second and third evolutionary lines from GenBank database. The A. ursinum plastome contains 90 predicted genes (81 unique) including 5 pseudogenes, while A. paradoxum has 88 predicted genes (79 unique) including 19 pseudogenes. The comparative analysis has revealed that the A. paradoxum plastome differs markedly from those of other species. Due to many deletions, the A. paradoxum plastome is the shortest of known for Allium species, being only 145,819 bp long. The most prominent distinctions are (1) a 4825 bp long local inversion that spans from the ndhE to the rpl32 gene in the small single copy region and (2) pseudogenization, or the loss of all NADH-genes. In contrast, the plastome of A. ursinum - a species from the first evolutionary line (as well as A. paradoxum) - resembles the Allium species of the second and third evolutionary lines, showing no large rearrangements or discrepancies in gene content. It is unclear yet whether only A. paradoxum was affected by some evolutionary events or its close relatives from both sect. Briseis and other sections of Amerallium were altered as well. We speculate the sunlight-intolerant, shade-loving nature of A. paradoxum and the impairment of the ndh genes in its plastome could be interrelated phenomena.

Erratum to "miR319, miR390, and miR393 Are Involved in Aluminum Response in Flax (Linum usitatissimum L.)".

[This corrects the article DOI: 10.1155/2017/4975146.].

miR319, miR390, and miR393 Are Involved in Aluminum Response in Flax (Linum usitatissimum L.).

Acid soils limit agricultural production worldwide. Major reason of crop losses in acid soils is the toxicity of aluminum (Al). In the present work, we investigated expression alterations of microRNAs in flax (Linum usitatissimum L.) plants under Al stress. Flax seedlings of resistant (TMP1919 and G1071/4_k) and sensitive (Lira and G1071/4_o) to Al cultivars and lines were exposed to AlCl3 solution for 4 and 24 hours. Twelve small RNA libraries were constructed and sequenced using Illumina platform. In total, 97 microRNAs from 18 conserved families were identified. miR319, miR390, and miR393 revealed expression alterations associated with Al treatment of flax plants. Moreover, for miR390 and miR393, the alterations were distinct in sensitive and resistant to Al genotypes. Expression level changes of miR319 and miR390 were confirmed using qPCR analysis. In flax, potential targets of miR319 are TCPs, miR390-TAS3 and GRF5, and miR393-AFB2-coding transcripts. TCPs, TAS3, GRF5, and AFB2 participate in regulation of plant growth and development. The involvement of miR319, miR390, and miR393 in response to Al stress in flax was shown here for the first time. We speculate that these microRNAs play an important role in Al response via regulation of growth processes in flax plants.

Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes.

BACKGROUND: The species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n = 30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes. RESULTS: High-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH. Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n = 30) could originate either as the result of hybridization of two diploid species (2n = 16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n = 16) and a diploid ancestor of modern L. narbonense (2n = 14). CONCLUSIONS: High-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.

Comparative analysis of inverted repeats of polypod fern (Polypodiales) plastomes reveals two hypervariable regions.

BACKGROUND: Ferns are large and underexplored group of vascular plants (~ 11 thousands species). The genomic data available by now include low coverage nuclear genomes sequences and partial sequences of mitochondrial genomes for six species and several plastid genomes. RESULTS: We characterized plastid genomes of three species of Dryopteris, which is one of the largest fern genera, using sequencing of chloroplast DNA enriched samples and performed comparative analysis with available plastomes of Polypodiales, the most species-rich group of ferns. We also sequenced the plastome of Adianthum hispidulum (Pteridaceae). Unexpectedly, we found high variability in the IR region, including duplication of rrn16 in D. blanfordii, complete loss of trnI-GAU in D. filix-mas, its pseudogenization due to the loss of an exon in D. blanfordii. Analysis of previously reported plastomes of Polypodiales demonstrated that Woodwardia unigemmata and Lepisorus clathratus have unusual insertions in the IR region. The sequence of these inserted regions has high similarity to several LSC fragments of ferns outside of Polypodiales and to spacer between tRNA-CGA and tRNA-TTT genes of mitochondrial genome of Asplenium nidus. We suggest that this reflects the ancient DNA transfer from mitochondrial to plastid genome occurred in a common ancestor of ferns. We determined the marked conservation of gene content and relative evolution rate of genes and intergenic spacers in the IRs of Polypodiales. Faster evolution of the four intergenic regions had been demonstrated (trnA- orf42, rrn16-rps12, rps7-psbA and ycf2-trnN). CONCLUSIONS: IRs of Polypodiales plastomes are dynamic, driven by such events as gene loss, duplication and putative lateral transfer from mitochondria.

The systematic position of Dryopteris blanfordii subsp. nigrosquamosa (Ching) Fraser-Jenkins within the genus Dryopteris Adans.

Dryopteris blanfordii (C.Hope) C.Chr. is a member of the Dryopteridaceae, growing in high altitude Picea or Abies forests (2900-3500 m) in China and India. Phylogenetic relationships between D. blanfordii subsp. nigrosquamosa and closely related species of Dryopteris were investigated using a combined analysis of multiple molecular data sets (the protein-coding region of rbcL and matK genes and intergenic spacers psbA-trnH, trnP-petG, rps4-trnS, trnL-trnF and rbcL-accD). An assumption about the position of D. blanfordii subsp. nigrosquamosa within Dryopteris was made by using the Maximum Likelihood and Bayesian Inference approach and chloroplast marker sequences of Dryopteris species from GenBank. The results demonstrated that Asian taxa D. blanfordii subsp. nigrosquamosa and D. laeta as well as two American species D. arguta and D. marginalis belong to the same clade, all four of them being part of Dryopteris section Dryopteris.

Protein Content of Circulating Nucleoprotein Complexes.

In the current study we have investigated the protein content of blood plasma deoxyribonucleoprotein complexes. The complexes were isolated using affinity chromatography with immobilized polyclonal anti-histone antibodies. Proteins were separated by SDS PAAGE and identified by MALDI-TOF mass-spectrometry. 111 and 56 proteins (excluding histones), respectively, were identified with a good score in deoxyribonucleoprotein complexes of healthy females and breast cancer patients. However, only four of these proteins were found in 30 % of all samples. Fourteen proteins previously described as tumor specific proteins were found in cancer patients whereas not one of them was found in healthy individuals. The data obtained demonstrate the involvement of different cellular and extracellular proteins in circulating cell-free DNA.

Identification, Expression Analysis, and Target Prediction of Flax Genotroph MicroRNAs Under Normal and Nutrient Stress Conditions.

Cultivated flax (Linum usitatissimum L.) is an important plant valuable for industry. Some flax lines can undergo heritable phenotypic and genotypic changes (LIS-1 insertion being the most common) in response to nutrient stress and are called plastic lines. Offspring of plastic lines, which stably inherit the changes, are called genotrophs. MicroRNAs (miRNAs) are involved in a crucial regulatory mechanism of gene expression. They have previously been assumed to take part in nutrient stress response and can, therefore, participate in genotroph formation. In the present study, we performed high-throughput sequencing of small RNAs (sRNAs) extracted from flax plants grown under normal, phosphate deficient and nutrient excess conditions to identify miRNAs and evaluate their expression. Our analysis revealed expression of 96 conserved miRNAs from 21 families in flax. Moreover, 475 novel potential miRNAs were identified for the first time, and their targets were predicted. However, none of the identified miRNAs were transcribed from LIS-1. Expression of seven miRNAs (miR168, miR169, miR395, miR398, miR399, miR408, and lus-miR-N1) with up- or down-regulation under nutrient stress (on the basis of high-throughput sequencing data) was evaluated on extended sampling using qPCR. Reference gene search identified ETIF3H and ETIF3E genes as most suitable for this purpose. Down-regulation of novel potential lus-miR-N1 and up-regulation of conserved miR399 were revealed under the phosphate deficient conditions. In addition, the negative correlation of expression of lus-miR-N1 and its predicted target, ubiquitin-activating enzyme E1 gene, as well as, miR399 and its predicted target, ubiquitin-conjugating enzyme E2 gene, was observed. Thus, in our study, miRNAs expressed in flax plastic lines and genotrophs were identified and their expression and expression of their targets was evaluated using high-throughput sequencing and qPCR for the first time. These data provide new insights into nutrient stress response regulation in plastic flax cultivars.

Gene expression profiling of flax (Linum usitatissimum L.) under edaphic stress.

BACKGROUND: Cultivated flax (Linum usitatissimum L.) is widely used for production of textile, food, chemical and pharmaceutical products. However, various stresses decrease flax production. Search for genes, which are involved in stress response, is necessary for breeding of adaptive cultivars. Imbalanced concentration of nutrient elements in soil decrease flax yields and also results in heritable changes in some flax lines. The appearance of Linum Insertion Sequence 1 (LIS-1) is the most studied modification. However, LIS-1 function is still unclear. RESULTS: High-throughput sequencing of transcriptome of flax plants grown under normal (N), phosphate deficient (P), and nutrient excess (NPK) conditions was carried out using Illumina platform. The assembly of transcriptome was performed, and a total of 34924, 33797, and 33698 unique transcripts for N, P, and NPK sequencing libraries were identified, respectively. We have not revealed any LIS-1 derived mRNA in our sequencing data. The analysis of high-throughput sequencing data allowed us to identify genes with potentially differential expression under imbalanced nutrition. For further investigation with qPCR, 15 genes were chosen and their expression levels were evaluated in the extended sampling of 31 flax plants. Significant expression alterations were revealed for genes encoding WRKY and JAZ protein families under P and NPK conditions. Moreover, the alterations of WRKY family genes differed depending on LIS-1 presence in flax plant genome. Besides, we revealed slight and LIS-1 independent mRNA level changes of KRP2 and ING1 genes, which are adjacent to LIS-1, under nutrition stress. CONCLUSIONS: Differentially expressed genes were identified in flax plants, which were grown under phosphate deficiency and excess nutrition, on the basis of high-throughput sequencing and qPCR data. We showed that WRKY and JAS gene families participate in flax response to imbalanced nutrient content in soil. Besides, we have not identified any mRNA, which could be derived from LIS-1, in our transcriptome sequencing data. Expression of LIS-1 flanking genes, ING1 and KRP2, was suggested not to be nutrient stress-induced. Obtained results provide new insights into edaphic stress response in flax and the role of LIS-1 in these process.

A clinical and molecular analysis of branchio-oculo-facial syndrome patients in Russia revealed new mutations in TFAP2A.

Branchio-oculo-facial syndrome (BOFS, OMIM# 113620) is a rare autosomal dominant disorder characterised by branchial cleft sinus defects, ocular anomalies and facial dysmorphisms, including lip or palate cleft or pseudocleft, and is associated with mutations in the TFAP2A gene. Here, we performed clinical analysis and mutation diagnostics in seven BOFS patients in Russia. The phenotypic presentation of BOFS observed in three patients showed high heterogeneity, including variation in its main clinical manifestations (linear loci of cervical cutaneous aplasia, ocular anomalies and orofacial cleft). In certain other cases, isolated ocular anomalies, or an orofacial cleft with accessory BOFS symptoms, were observed. In five BOFS patients, conductive hearing loss was diagnosed. Direct sequencing of the coding region of the TFAP2A gene revealed missense mutations in four BOFS patients. One patient was observed to have a previously described mutation (p.Arg251Gly), while three patients from two families were found to have novel mutations: p.Arg213Ser and p.Val210Asp. These novel mutations were not present in healthy members of the same family and therefore should be classified as de novo.

Excess fertilizer responsive miRNAs revealed in Linum usitatissimum L.

Effective fertilizer application is necessary to increase crop yields and reduce risk of plant overdosing. It is known that expression level of microRNAs (miRNAs) alters in plants under different nutrient concentrations in soil. The aim of our study was to identify and characterize miRNAs with expression alterations under excessive fertilizer in agriculturally important crop - flax (Linum usitatissimum L.). We have sequenced small RNAs in flax grown under normal and excessive fertilizer using Illumina GAIIx. Over 14 million raw reads was obtained for two small RNA libraries. 84 conserved miRNAs from 20 families were identified. Differential expression was revealed for several flax miRNAs under excessive fertilizer according to high-throughput sequencing data. For 6 miRNA families (miR395, miR169, miR408, miR399, miR398 and miR168) expression level alterations were evaluated on the extended sampling using qPCR. Statistically significant up-regulation was revealed for miR395 under excessive fertilizer. It is known that target genes of miR395 are involved in sulfate uptake and assimilation. However, according to our data alterations of the expression level of miR395 could be associated not only with excess sulfur application, but also with redundancy of other macro- and micronutrients. Furthermore expression level was evaluated for miRNAs and their predicted targets. The negative correlation between miR399 expression and expression of its predicted target ubiquitin-conjugating enzyme E2 gene was shown in flax for the first time. So we suggested miR399 involvement in phosphate regulation in L. usitatissimum. Revealed in our study expression alterations contribute to miRNA role in flax response to excessive fertilizer.

Retrotransposon-based molecular markers for analysis of genetic diversity within the Genus Linum.

SSAP method was used to study the genetic diversity of 22 Linum species from sections Linum, Adenolinum, Dasylinum, Stellerolinum, and 46 flax cultivars. All the studied flax varieties were distinguished using SSAP for retrotransposons FL9 and FL11. Thus, the validity of SSAP method was demonstrated for flax marking, identification of accessions in genebank collections, and control during propagation of flax varieties. Polymorphism of Fl1a, Fl1b, and Cassandra insertions were very low in flax varieties, but these retrotransposons were successfully used for the investigation of Linum species. Species clusterization based on SSAP markers was in concordance with their taxonomic division into sections Dasylinum, Stellerolinum, Adenolinum, and Linum. All species of sect. Adenolinum clustered apart from species of sect. Linum. The data confirmed the accuracy of the separation in these sections. Members of section Linum are not as closely related as members of other sections, so taxonomic revision of this section is desirable. L. usitatissimum accessions genetically distant from modern flax cultivars were revealed in our work. These accessions are of utmost interest for flax breeding and introduction of new useful traits into flax cultivars. The chromosome localization of Cassandra retrotransposon in Linum species was determined.

Human gut microbiota community structures in urban and rural populations in Russia.

The microbial community of the human gut has a crucial role in sustaining host homeostasis. High-throughput DNA sequencing has delineated the structural and functional configurations of gut metagenomes in world populations. The microbiota of the Russian population is of particular interest to researchers, because Russia encompasses a uniquely wide range of environmental conditions and ethnogeographical cohorts. Here we conduct a shotgun metagenomic analysis of gut microbiota samples from 96 healthy Russian adult subjects, which reveals novel microbial community structures. The communities from several rural regions display similarities within each region and are dominated by the bacterial taxa associated with the healthy gut. Functional analysis shows that the metabolic pathways exhibiting differential abundance in the novel types are primarily associated with the trade-off between the Bacteroidetes and Firmicutes phyla. The specific signatures of the Russian gut microbiota are likely linked to the host diet, cultural habits and socioeconomic status.

MALINA: a web service for visual analytics of human gut microbiota whole-genome metagenomic reads.

MALINA is a web service for bioinformatic analysis of whole-genome metagenomic data obtained from human gut microbiota sequencing. As input data, it accepts metagenomic reads of various sequencing technologies, including long reads (such as Sanger and 454 sequencing) and next-generation (including SOLiD and Illumina). It is the first metagenomic web service that is capable of processing SOLiD color-space reads, to authors' knowledge. The web service allows phylogenetic and functional profiling of metagenomic samples using coverage depth resulting from the alignment of the reads to the catalogue of reference sequences which are built into the pipeline and contain prevalent microbial genomes and genes of human gut microbiota. The obtained metagenomic composition vectors are processed by the statistical analysis and visualization module containing methods for clustering, dimension reduction and group comparison. Additionally, the MALINA database includes vectors of bacterial and functional composition for human gut microbiota samples from a large number of existing studies allowing their comparative analysis together with user samples, namely datasets from Russian Metagenome project, MetaHIT and Human Microbiome Project (downloaded from http://hmpdacc.org). MALINA is made freely available on the web at http://malina.metagenome.ru. The website is implemented in JavaScript (using Ext JS), Microsoft .NET Framework, MS SQL, Python, with all major browsers supported.

Impact of recombination on polymorphism of genes encoding Kunitz-type protease inhibitors in the genus Solanum.

BACKGROUND: The group of Kunitz-type protease inhibitors (KPI) from potato is encoded by a polymorphic family of multiple allelic and non-allelic genes. The previous explanations of the KPI variability were based on the hypothesis of random mutagenesis as a key factor of KPI polymorphism. RESULTS: KPI-A genes from the genomes of Solanum tuberosum cv. Istrinskii and the wild species Solanum palustre were amplified by PCR with subsequent cloning in plasmids. True KPI sequences were derived from comparison of the cloned copies. "Hot spots" of recombination in KPI genes were independently identified by DnaSP 4.0 and TOPALi v2.5 software. The KPI-A sequence from potato cv. Istrinskii was found to be 100% identical to the gene from Solanum nigrum. This fact illustrates a high degree of similarity of KPI genes in the genus Solanum. Pairwise comparison of KPI A and B genes unambiguously showed a non-uniform extent of polymorphism at different nt positions. Moreover, the occurrence of substitutions was not random along the strand. Taken together, these facts contradict the traditional hypothesis of random mutagenesis as a principal source of KPI gene polymorphism. The experimentally found mosaic structure of KPI genes in both plants studied is consistent with the hypothesis suggesting recombination of ancestral genes. The same mechanism was proposed earlier for other resistance-conferring genes in the nightshade family (Solanaceae). Based on the data obtained, we searched for potential motifs of site-specific binding with plant DNA recombinases. During this work, we analyzed the sequencing data reported by the Potato Genome Sequencing Consortium (PGSC), 2011 and found considerable inconsistence of their data concerning the number, location, and orientation of KPI genes of groups A and B. CONCLUSIONS: The key role of recombination rather than random point mutagenesis in KPI polymorphism was demonstrated for the first time.